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The faint outskirts of the blue compact galaxy Haro 11: is 
there a red excess? 
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ABSTRACT 

Previous studies of the low surface brightness host of the blue compact galaxy (BCG) 
Haro 11 have suggested an abnormally red color of V— K = 4.2±0.8 for the host galaxy. 
This color is inconsistent with any normal stellar population over a wide range of stellar 
metallicities (Z — 0.001-0.02). Similar though less extreme host colors have been 
measured for other BCGs and may be reconciled with population synthesis models, 
provided that the stellar metallicity of the host is higher than that of the ionized gas 
in the central starburst. We present the deepest V and K band observations to date 
of Haro 11 and derive a new V — K color for the host galaxy. Our new data suggest 
a far less extreme colour oi V — K = 2.3 ± 0.2, which is perfectly consistent with 
the expectations for an old host galaxy with the same metallicty as that derived from 
nebular emission lines in the star-forming center. 
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1 INTRODUCTION 

BCGs are gas-rich low-luminosity galaxies of particular in- 
terest since many of them have very hig h star formation 
rates (SFR) and low chemical abundances (|Kunth fc OstlirJ 
2000). They are thus reminiscent of high-redshift star- 
bursting young galaxies. This, together with their close 
proximity to us, makes BCGs useful test objects, suit- 
able for gaining insight into galaxy formation and intense 
star formation (SF). Aside from the bright central star- 
burst, deep photometric optical and near- infrared (NIR 
st udies (e.g. lLoose fc ThuarJ Il986l: iDoublier et all m^T 



Papaderos et al.lll99rj|cair6s et alj|200ll; IBergvall fc OstlirJ 
20021 : lAmorfn et alj|200l 120091 : iNoeske et al.ll2005l ) have re- 
vealed the presence of another component of low surface 
brightness (LSB). The LSB component hosts the central 
starburst and is therefore referred to as a LSB host galaxy. 
The host is only visible close to the outskirts of BCGs since 
at small radial distances the intensity from the central star- 
burst outshines the contribution from the LSB component. 
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Many photometric studies have concentrated on reaching 
faint levels in an attempt to isolate the LSB component and 
characterize it. Some have suggested that the LSB compo- 
nent has unusual propertie s such as an extreme red excess 
in the optical/N IR colors. IBergvall fc OstlirJ |2002l . here- 
after BO02) and IBergvall et al.l |2005t ) present a sample of 
10 BCGs with progressively redder colors towards t he out- 
skirts of the host. BO02 and fZackriss orTet al.l (|2006l . here- 
after Z06) argue that these colors are much too red to be 
reconciled with a normal stellar population. 

Extended faint structures of very red colors have 
also been detected in the outskirts of disk galaxies 
and, similarly to BCGs, these structures display colors 
much too red to be due to stellar populations with 
reasonable metallicities and normal initial mass fu nc- 
tions (IMFs). In 2004, IZibetti. White fc BrinkmanrJ re- 
ported a faint structure with a notable red excess 
around the stacked disk of 1047 edge-on spirals from 
the Sloan Digital Sky Survey (SDSS). In the same 
year, IZibetti fc Ferguson! reported a red structure around 
a disk galaxy in the Hubble Ultra Deep Field (HUDF). Sim- 
ilarly, IBergvall. Zackrisson fc Caldwell! |2009l ) find a red ex- 
cess in the outskirts of disk galaxies by stacking 1510 edge-on 
LSB galaxies from the SDSS. 

These detections have been collectively dubbed "red halos" 
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and a number of possible explanations have been proposed. 
Some detections have been f alsified and ascr ibed to instru- 
mental effects. For instance, Ide Jongl (|2008h demonstrated 
that the HUDF red halo is likely due to reddening by the 
far wings of the point spread function (PSF). Most detec- 
tions, however, have persisted and even given rise to exotic 
scenarios. BO02 proposed that the extreme colors seen in 
their sample could be explained by a metal rich (> Zq) 
halo stellar population, in contrast to the low nebular chem- 
ical abundances. Z06 argued that the only stellar popula- 
tion which explains the published colors of the halos of both 
disk galaxies and BCGs is one with an extremely bottom- 
heavy (dN/dM oc M~ a with a = 4.50) IMF. If proven 
to be a common component of galaxies of various Hubble 
types, this would imply that faint low-mass stars are so nu- 
merous in these galaxies that they may significantly con- 
tribute to the baryons still missing from inventories in the 
low-redshift Universe. Another effect which may produce red 
halos around both disk galaxies and BCGs is that the surface 
brightness profiles used to derive the colors have been over- 
subtracted due to failure to ac count for extinction by the ex- 
tragalactic background light jZackrisson. Micheva fc Ostlinl 
l2009h . 

If one is willing to trade off a common explanation for all 
red halo detections with individual explanations applicable 
to certain galaxies and not to others then mo st of the red 
halos of the BO02 and iBergvall et"afl (120051 ) BCGs sam- 
ple can be modeled with a Salpeter IMF, with the caveat 
that the required stellar metallicity is in some cases much 
higher th an what is observed fr om the gas in the central 
starburst dZackrisson et all 120091 '). Among the sample, one 
of the biggest, most massive and most luminous BCGs is 
Haro 11, with a very high SFR of 18-20 M yr _1 and the 
reddest color of V — K = 4.2 ± 0.8, which cannot be recon- 
ciled with any normal stellar population. 
Throughout this paper we assume a Hubble parameter of 
Ho = 75 km s" 1 Mpc -1 and a distance to Haro 11 of 
82 Mpc, with a redshift of z = 0.020598. In §[2] we present 
the deepest observations to date for Haro 11 in the V and 
K passbands and a new and more accurate measurement 
of the V — K color. We are pushing the limits of surface 
photometry down to isophotes 10000 times fainter than the 
brightness of the night sky in K and we therefore examine 
the reliability of the obtained profiles in § 12.31 and § 12.41 
Our results for the V — K color of Haro 11 are presented in 
§ Discussion and summary can be found in § [3] and § 
respectively. 



2 DATA 

Table [T] summarizes the available data used in this paper 
and co mpares with the data published bv lBergvall fc Ostlinl 
(2002). The V band data was obtained in 2008 at the 
Nordic Optical Telescope (NOT) with the MOSaic CAmera 
(MOSCA), which has a scale of 0.217 arcsec/pixel. The K 
band data was obtained in 2005 at the ESO NTT 3.58m tele- 
scope with Son OF Isaac (SOFI), which has a scale of 0.288 
arcsec/pixel. The raw frame size in the optical (MOSCA) is 
7.7 arcmin, in the NIR (SOFI) 4.9 arcmin. The NIR images 
are dithered with large offsets. Due to the large field of view 
of the instrument and the comparatively much smaller tar- 



get size, this is the optimal observational strategy for this 
object. 

2.1 Reductions 

The data were reduced in the following way. The raw images 
in both bands were cleaned from bad pixels, bias subtracted 
in the V band, pair subtracted in the K band using the 
sky median as a scaling factor, trimmed, flatfielded with an 
illumination-corrected normalized masterflat, and sky sub- 
tracted with a fiat surface interpolated with a first order 
polynomial. In the V band the final stacked image was cali- 
brated with secondary standard stars against the published 
data by B002. In the K band each sky-subtracted frame 
was calibrated against 2MASS stars located in the frame, 
whereby a minimum of five 2MASS stars per frame were 
used to obtain its zeropoint. The calibrated K band data 
was therefore independent of the photometric conditions, 
even though an examination of the temporal zeropoint vari- 
ation throughout the observations revealed a scatter of only 
~ 0.049 magnitudes. In both bands, images were registered 
using Pyraf GEOMAP/GEOTRAN. When shifting and ro- 
tating data we interpolated with a 5th order polynomial in 
order to avoid the Moire effect caused by the default bilinear 
interpolation. The images were median and average com- 
bined without the use of a rejection algorithm. Comparison 
between the median and average stacked images revealed an 
insignificant difference in the aperture photometry of point 
sources on the order of ~ 0.02 magnitudes. Median stacking 
has the added advantage of filtering out any lingering cosmic 
rays and bad pixels, therefore the analysis in this paper was 
carried out on the median stacked images in both bands. Sky 
subtraction was performed on the final stacked image in each 
band in order to remove any residual sky. When perform- 
ing sky subtraction we neglected the effects of dust extinc- 
tion of the extrag alactic background light in the out skirts 
of the host galaxy |Zackrisson. Micheva fc 6stlinll2009h and 
assumed that the sky over the target can be interpolated 
from sky regions around the target. The final reduced im- 
ages in V and K have an effective field of view smaller than 
the field of view of the detectors in both bands, which is a 
consequence of the median stacking of dithered images. 

2.2 Sky estimation 

The shape of a radial color profile is sensitive to the slope 
of the two surface brightness profiles from which it has been 
derived. In turn, the slope of the surface brightness profiles 
is sensitive to any systematic errors in the sky subtraction, 
especially at large radial distances and very faint isophotal 
levels. In particular, over- or undersubtracting the sky back- 
ground would lead to spurious color gradients. It is there- 
fore important to have control over all aspects of the sky 
background estimation. We briefly outline the major steps 
involved in this procedure. 

For each frame the sky is approximated by fitting a flat sur- 
face to re gions free from so urces ("sky regions") with the 
PyMidas (|Hook et all 120061 ) procedure FitFlat. The place- 
ment and size of these sky regions are very important - one 
aims to have suitably large sky regions uniformly distributed 
over the frame in order for the interpolated sky surface to 
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Old data New data 

V K' V Ks 

Year 1984 1993 2008 2005 

Exposure Time 15 min 66 min 40 min 60 min 

Instrument ESO 2.2m ESO 2.2m IRAC2 NOT 2.56m MOSCA ESO NTT 3.58m SOFI 

///pixel 0.362 0.475 0.217 0.288 

FoV(//) 178x178 136x136 462x462 295x295 



Table 1. Observational data summary for Haro 11 (ESO350-IG038). The effective exposure time for new 
V (new K) data is 3.3 (2.4) times longer than for the old V (old K) data. Haro 11 has an apparent angular 
diameter of ~ 1 arcmin out to a surface brightness of fiy = 26.5 mag arcsec -2 . All observations are in 
dithering mode. 



be a well-sampled approximation of the sky. The size of the 
sky regions together with the order of the interpolating func- 
tion determine whether the interpolated surface will mimic 
the large scale or the small scale structure of the sky. It 
is difficult to apply a reasonable model to small scale sky 
structures. We have chosen to remove only the large scale 
structures, but include the remaining small scale sky fluctu- 
ations in the error analysis (§ 12. 3p . 

In order to model the large scale sky structure of a frame, 
the area of each sky region is not allowed to drop below 
600 pixels 2 and the interpolating function is a polynomial 
of first order in both x and y. The position and size of the 
sky regions are determined automatically by the following 
algorithm. All positive and negative sources are masked out, 
giving a binary mask for each image. Initially, a very coarse 
fixed grid is placed over the mask image. Defining sky regions 
over a frame thus translates to iterating over gridboxes and, 
at each iteration, determining the maximum size of the sky 
region that can be placed inside each gridbox. The size of 
the sky region depends on the presence of sources. If sources 
are absent, the sky region takes the shape and size of the 
entire gridbox. If sources are present, a box of smaller size 
is generated and tested for sources. Smaller boxes are first 
generated by splitting the gridbox at an iteratively moving 
point along one axis while keeping the size of the box along 
the other axis at maximum. These boxes are thus always an- 
chored at two grid vertices along the same axis. If all of the 
resulting boxes test positive for sources new boxes are gen- 
erated by moving a strip of fixed width along one axis and 
maximum height in the other axis. These boxes may or may 
not be anchored to any vertices but always touch two oppos- 
ing gridbox axes. If this too fails, new boxes are generated 
by varying their size along both axes. These boxes are not 
adjacent to vertices or to gridbox axes and instead "float" 
in the gridbox interior. If still no boxes free from sources are 
found the gridbox is simply split in four quadrants along the 
middle of both its axes, resulting in four boxes of a quar- 
ter size. If these quarter boxes also test positive for sources, 
then no sky region can be defined over this gridbox and the 
algorithm continues to the next gridbox. 

After going through the entire grid each quadrant of the 
frame is checked for the number of obtained boxes. If the 
number of sky boxes in any quadrant is less than 20, the 
grid is refined and iteration starts over. The algorithm is 
such that the exact position, shape and size of the present 
sources is irrelevant. The imposed limit on the minimum 
number of sky boxes in each frame quadrant ensures that 
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Figure 3. Radial error composition. Uncertainties in mag- 
nitudes vs. distance in arcseconds from the center of Haro 11 in 
the V band (upper panel) and the K band (lower panel). The 
uncertainty in the sky, o^kyi is marked with filled squares; the 
standard deviation of the mean, ctsdoMj with open circles; the 
total error, <r tot , with an uninterrupted line. The large contribu- 
tion of ctsdom at small radii is due to resolved physical structure 
present at the center of Haro 11. For both bands c sky dominates 
the total error at large radii. 

the tiepoints for the interpolated surface are well distributed 
to cover the entire frame. 

We have also carried out tests to evaluate the performance 
of the automated sky box placement algorithm. These con- 
sisted of varying the box shape, box size, and the number 
of boxes per quadrant, and comparing the effect that would 
have on the PyMidas FitFlat sky estimation. We conclude 
that the algorithm is stable and reliable in all of our test 
cases with the caveat that spatially non-linear low frequency 
sky variations are intrinsically impossible to perfectly model 
with a first order polynomial regardless of the choice in 
placement and size of sky boxes. While we have used PyMi- 
das FitFlat for the sky subtraction itself, a test with IRAF 
IMSURFIT gives nearly identical results, provided that the 
same sky box distribution is used in both cases. 

2.3 Surface brightness profiles 

The surface brightness profiles in the right panel of Fig. [5] 
are obtained by integrating in elliptical rings starting from 
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Figure 1. New V and K band data for Haro 11. Left panel: NOT MOSCA V band with scale 0.217"/plxel, 40 minute exposure, 1 
chip. The image size is 2.7' X 2.3'. Right panel: ESO NTT SOFI K band with scale 0.288" /pixel, 1 hour exposure. The image size is 
2.3' X 2.2'. For display purposes the K band image (right) has been scaled to match the size of the V band image (Note that the images 
have different center coodinates). Bad pixels, sources and cosmic rays in both images are masked out. The larger size of the images 
significantly increases the area available for the estimation of the sky background. North is up, east is to the left. The large elliptical 
mask in the V band image is due to a bright blue star. 




Figure 2. Comparison of the K band data for Haro 11. Left panel: the data used in the analysis of iBergvall fc OstlirJ l|2002l l. The 
image size is ~l'xl' with a scale of 0.475 arcsec/pixel. Right panel: the central ~l'xl' of the new K band data with a scale of 0.288 
arcsec/pixel. North is up, east is to the left. The new image is of better quality, with higher resolution and a flatter background. 



the center of the galaxy. The faint outskirts of Haro 11 have 
an almost circular shape in V, however the central starburst 
is very bright and it is unclear at what radii its contribution 
to the overall intensity becomes negligible. The shape and 
center of the underlying host is therefore better estimated 



from the K band, where the contribution from the starburst 
quickly drops off with increasing radius. The coordinates 
of the center and the parameters of the ellipse are chosen 
by fitting isophotes to the K band image with IRAF EL- 
LIPSE. For robustness, different centering positions, inclina- 
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Table 2. Error composition of the host defined over the radial 
range 6-10 kpc. All error values are in magnitudes. cr a ^ y is the 
uncertainty in the sky background, estimated from the standard 
deviation of the mean values of 275 sky regions in K and 280 
sky regions in V. The range of sizes of the sky regions is (12 X 
13, 16x21, 15x27, 21x21, 27x27) pixels, ctsdom is the standard 
deviation of the mean of the flux inside the elliptic ring between 
6 and 10 kpc. otot is the composite error, with <T s ky and o"SDOM 
added in quadrature. 



tions and position angles (P.A.) have been tested. There is 
no discernible difference in the radial color profiles obtained 
from integration with ellipticity and P.A. parameters in the 
range e G (0.20,0.34), P.A. £ (103°, 120°), the latter mea- 
sured from North through East. The range of these values 
covers all elliptical isophotes returned by IRAF ELLIPSE 
at surface brightness levels of fiK = 21-23.5 mag arcsec -2 , 
where the signal from the host galaxy dominates the lumi- 
nosity output. These values are similar to the ones used by 
BO02 (e = 0.20, P.A. = 120°). Once integration starts, the 
parameters of the elliptic rings and the step size along the 
major axis are kept constant. The profiles are sampled with 
a step size of 5 pixels. 

At faint isophotal levels it is the uncertainty in the zero 
level of the sky background, a s k y , that largely dominates 
the error (Fig. [3j . This uncertainty is estimated by masking 
out all sources, measuring the mean intensity inside square 
apertures, and calculating the standard deviation of these 
means. The a 3 k y obtained in this fashion is by definition 
always numerically greater than the standard deviation of 
the mean sky level and thus represents a conservative esti- 
mate of the sky uncertainty. The apertures are of the same 
order in size as the area of the smallest elliptic ring inside 
of the LSB host, measured at /ik ~ 21.5 mag arcsec -2 . An- 
other source of uncertainty included in the error analysis is 
the uncertainty in the mean flux level of each ring, which 
is a composite error of the poisson noise and the intrinsic 
flux scatter across each ring, given by the standard devia- 
tion of the mean (SDOM) . Instead of isolating the poisson 
noise in each ring, we use the SDOM in the rings as an upper 
limit for the poisson noise, because the reduction process in- 
variably destroys the initial poisson distribution of the raw 
frames. To obtain a more accurate value, one would have 
to obtain the poisson noise per ring on the raw frames and 
then propagate those values throughout the reduction pro- 
cess. For the purposes of our analysis, however, the SDOM 
of each ring on the final reduced frames is an acceptable 
upper limit. Thus, the combination of the sky uncertainty 
and the SDOM gives a conservative representation for the 
uncertainties in our results, summarized in Table [2] 
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Figure 4. Surface brightness profiles of synthetic exponential 
disks. Filled circles: steep disk with scale length h T = 2.1 kpc; 
open circles: Haro 11-like disk with h T = 2.7 kpc; filled squares: 
flat disk with h T = 4.0 kpc. Circles and squares are used for 
the profiles of the synthetic disks recovered after reduction and 
calibration. The dashes lines are the analytical profiles of each 
synthetic disk on an empty frame with no noise, /ig = 23 
mag arcsec -2 is the faintest isophote that can be reached by all 
three exponential disk profiles with negligible random scatter and 
no systematic effects, see text for justification. This is marked by 
the gray-shaded area in the figure. 



2.4 Setting limits on systematic errors in surface 
brightness profiles 

At low signal-to-noise levels the interpretation of a surface 
brightness profile is subject to some ambiguity. Observed 
profile trends could either be real or due to the propaga- 
tion of systematic errors that have not been accounted for 
in the reduction process. It is also not always obvious at 
what surface brightness levels the sky noise starts domi- 
nating the behavior of the profile. When pushing surface 
photometry down to very faint isophotes it is necessary to 
identify the maximum radial range within which the data 
are not dominated by any systematic errors that could have 
been introduced by the reduction process or by the random 
fluctuations in the sky background. Attempting to ascertain 
the effects of reduction on the data by analyzing the sur- 
face brightness profile of a real physical target is not always 
practical since the true shape and slope of the profile is often 
not known in advance. 

A better way to perform such a test is to insert a synthetic 
galaxy with an a priori known surface brightness profile in 
all raw frames, then reduce and calibrate with the same 
pipeline that was used on the real data. The easiest profile 
to identify is that of an exponential disk, since it appears 
as a straight line in a log-lin plot. Analyzing the surface 
brightness profile of a reduced exponential disk will reveal 
the presence of any systematic errors, as well as give an 
indication at what surface brightness levels the sky noise 
starts dominating the profile. The surface brightness level, 
at which the reduced synthetic disk profile starts deviating 
from the expected straight line in a log-lin plot marks the 
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faintest isophote that can be reached with the quality of the 
data and the accuracy of the pipeline. 

It is the K band, rather than the V band, that has the nois- 
iest sky and the largest uncertainty in the sky level, so a 
V — K profile will inevitably be subject to limits imposed 
solely by the depth and quality of the K band data. We 
therefore perform the exponential disk test only on the K- 
band data. For robustness the test is carried out three times, 
each with different positions and scale lengths of the disks. 
Fig. [4] shows a plot of the surface brightness profiles of the 
three reduced exponential disks. The scale lengths are cho- 
sen to be steeper, approximately equal to, and flatter than 
the scale length of the Haro 11 LSB host. The latter is ob- 
tained by a least-squares fit to the Haro 11 profile in a range 
of jj,K = 21.5-23.0 mag arcsec -2 , which is well away from the 
central starburst. This gives a Haro 11 scale length of w 2.6 
kpc and a y- intercept of /io,A' = 19.6 mag arcsec -2 . The 
central surface brightness is kept constant for all three syn- 
thetic disks at (iq.k ~ 19.5 mag arcsec -2 . Scale lengths of 
2.1, 2.7 and 4.0 kpc are chosen for the steeper, the approx- 
imately equal, and the flatter synthetic disks respectively. 
The shortest and longest scale lengths are chosen to cor- 
respond to radial surface brightness profiles that are ~ l m 
fainter, respectively ~ l m brighter than a profile with a scale 
length of 2.7 kpc at a distance of r = 25 arcsec where the 
surface brightness is [ik = 23 mag arcsec -2 . This allows 
the results of the test to be applicable to a range of profile 
slopes, without the need for great accuracy in the determi- 
nation of the true scale length of the host galaxy. Indeed, 
measurements of the Haro 1 1 scale length from various loca- 
tions and lengths within /ik = 20.5-23.8 mag arcsec -2 give a 
slightly different value, however the difference is never found 
to be large enough to displace the profile from within the 
envelope created by the steep and flat exponential disks in- 
cluded in the test. The K band surface brightness profile for 
Haro 11 would have to be off by more than ±l m than what 
we observe at the radial distance of r = 25 arcsec, in order 
for this test to significantly over- or underestimate the K 
band profile depth that can be reached. 
Along with the surface brightness profiles of the three re- 
duced exponential disks, Fig. [4] also shows the analytical 
profiles of each disk on an empty frame. In the absence 
of any systematic and random uncertainties or in regions 
where these are insignificant, the reduced exponential disk 
profiles will exactly follow their analytical counterparts. The 
steeper and the Haro 11-like disk profiles are consistent with 
their respective analytical profiles out to r ~ 20.7 arcsec 
and r ~ 27.6 arcsec, respectively, both reaching a depth 
of hk ~ 23.5 mag arcsec -2 with a maximum deviation 
from their analytical profiles of 5 < 0.15 magnitudes at 
these radii. The flat disk, however, deviates from its analyt- 
ical counterpart by more than 5 « 0.3 magnitudes beyond 
fAK ~ 23.0 mag arcsec -2 , r > 32.5 arcsec. We note that 
accepting a deviation of 5 < 0.3 magnitudes would make 
this profile consistent with its analytical counterpart down 
to fiK ~ 24.5 mag arcsec -2 , r ~ 46.5 arcsec, however the 
corresponding maximum deviations of the steep and Haro 
11-like profiles are then S < 0.45 and S < 0.75 magnitudes, 
respectively, which we deem to be unacceptably large. 
Analyzing the exact behavior of different exponential disks 
at extremely low surface brightness levels is, however, be- 
yond the scope of this paper. We are therefore content with 



using the results of this test to set a conservative limit on 
the depth of the K band data and we conclude that a sur- 
face brightness level of [lk = 23 mag arcsec -2 can clearly be 
reached for all three profiles with negligible random scatter 
(8 < 0.15) and no discernible systematic effects. We adopt 
the limit of )j,k = 23 mag arcsec -2 as the faintest isophote 
that can be reached in the K band with the quality of the 
Haro 11 data and with the current pipeline. 
Using a synthetic disk with a Sersic index n > 1 will not 
have an effect on the outcome or the interpretation of this 
test. At the relevant radial distances away from the center 
all profiles with reasonable values for n are essentially par- 
allel to the exponential disk. Moreover, the exact shape of 
the analytical profile is irrelevant. It is only important at 
what radial distance from the center the profile shape can 
no longer be recovered with acceptable scatter and no sys- 
tematic effects. 



3 RESULTS 

The V and K surface brightness profiles are presented in 
Fig. E] together with the profiles published by BO02. The V 
profile is well-behaved and with reasonably small errorbars 
down to nv ~ 28.5 mag arcsec -2 . For the K band, however, 
the synthetic disk tests in § !2.4l suggest that a K band profile 
obtained from the current data with the current pipeline is 
reliable down to fj,K ~ 23 mag arcsec -2 . At fainter isophotes 
the profile is likely to be dominated by systematic errors 
introduced by the reduction procedure. For the Haro 11 K 
band profile in Fig. \5\ (ik = 23 mag arcsec -2 is reached at 
a radial distance of r ~ 10 kpc, or ~ 25 arcsec, from the 
center of the galaxy. This radius constitutes the maximum 
radius that we consider in the consecutive calculation of the 
total color for the Haro 11 LSB host galaxy. 
Fig. [5] shows the V — K radial color profile together with the 
profile published by BO02. Both profiles are uncorrected for 
dust extinction and can be directly compared. The two pro- 
files display different behavior. The BO02 profile shows no 
indication of breaking the trend of becoming progressively 
redder outwards of r 5 arcsec, reaching V — K ~ 3 at 
r ~ 15 arcsec, while the profile of the new and deeper data 
reddens from r ~ 5 arcsec to r ~ 15 arcsec only to flatten 
out and stabilize at V — K ~ 2.3 outwards of r ^ 15 arcsec. 
For the calculation of the total V — K halo color we use 
the radial range r ~ 5-10 kpc. The K band profile changes 
slope at r ~ 5 kpc and seems consistent with a constant 
slope out to r ~ 13 kpc. This suggests that the contribution 
from the central starburst in the K band is negligible already 
at r ~ 5 kpc and beyond this radius the profile is instead 
dominated by the LSB host. We therefore take ro — 5 kpc as 
the starting radius for the calculation of the total color. For 
the maximum radius we take the limit of r max — 10 kpc, 
as previously justified in § 12.41 Over this radial range, we 
obtain a total color of V — K — 2.3 ± 0.2. This is 2.3<r away 
from the BO02 value of V - K = 4.2 ± 0.8. 
It should be noted that the total V — K color reported in 
BO02, V — K = 4.2 ± 0.8, is measured out to a distance 
of r ~ 12.3 kpc. This is outside of the maximum radius 
limit we impose on our data, r max ~ 10 kpc. For direct 
comparison of the two total color measurements, we also 
measure the total color out to a radius of r ~ 12.3 kpc as 
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Figure 5. Surface brightness profiles of Haro 11 in V (triangles) and K (circles). Left panel: profiles by BO02 with e = 0.20, P. A. = 120°. 
The error bars represent the uncertainty in the sky level estimated to be the spread in mean values of a sample of "empty" sky regions 
uniformly distributed over the frames.; Right panel: profiles derived in this work with e = 0.32, P. A. = 103°. The error bars represent 
the uncertainty in the sky level, (T s k y , together with the uncertainty of the mean flux value of each radius, ctsdom- The bumpy feature 
inside r < 5" is due to Haro ll's 3 brightest knots l|Cumming et al.ll2009h . neither of which is at the center of profile integration. The 
horizontal line at = 23 mag arcsec - 2 marks the faintest isophote beyond which the K band profile of an exponential disk systematically 
deviates from a straight line by more than 0.3 magnitudes. For comparison, the same line is marked in the left panel. 
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Figure 6. V — K color profiles of Haro 11. Left panel: profile by BO02 with a step size of 0.95". The error bars represent the uncertainty 
in the zero level of the sky. No correction for dust reddening is applied. The profile becomes progressively redder with increased radius.; 
Right panel: profile derived in this work with a step size 1.44". The error bars represent the uncertainty in the sky level, (X^y, together 
with the uncertainty of the mean flux value of each radius, ctsdom- No correction for dust reddening is applied. The bumpy feature inside 
r < 5" is due to Haro ll's 3 brightest knots, neither of which is at the center of profile integration. At a radius of r ~ 14"the color 
stabilizes at V — K ~ 2.3. The gray-shaded area corresponds to surface brightness magnitudes fainter than fj,fc = 23 mag arcsec -2 and 
is not included in the calculation of the total V — K color. For comparison, the same range is marked in the left panel. 



in BO02, and obtain the slightly higher value for the total 
color of V — K — 2.4±0.5, which is 2a away from the total 
color of BO02. We cannot justify the use of data points 
out to such large radii in the total V — K color because 
our exponential disk test with a Haro 11-like disk shows 
a systematic brightening of the K band surface brightness 
profile beyond r > 11 kpc (Fig. [4|. In the absence of a 
similar systematic effect in the V band at the same radii, the 



brightening of the K band profile will result in a spurious 
reddening of the V — K radial profile. 

We further note that a simple least-squares fit of the Haro 
11 profile in the range of [ik = 21.5-23.0 mag arcsec -2 gives 
a scale length of the host population of ~ 2.6 kpc in the 
K band with an extrapolated central surface brightness of 
(io, k — 19.6 mag arcsec -2 . 
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Figure 7. Marigo ct al. (2008) stellar evolutionary tracks for dif- 
ferent met allicities with halo co lor data for all 10 BCGs of the 
BO02 and lBergvall et alj j2005l l sample. The old color V - K = 
4.2 ± 0.8 for Haro 11 is marked with a filled square. The tracks 
clearly cannot reproduce this data point. 
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Figure 8. Marigo ct al. ( 200^) stellar evolutionary tracks for Z = 
0.001 and Z = 0.004 with the new Haro 11 halo color of V-K = 
2.3±0.2 marked with a filled triangle. The data point is consistent 
with an old stellar population of low metallicity, Z = 0.001. 



4 DISCUSSION 

In this section we examine how our results compare with 
previously published measurements and comment on the im- 
plications of the new Haro 11 host color for the underlying 
stellar population. 



tion of low metallicity. This has lead to somewhat con- 
troversial sc enarios for the nature of the underlying stellar 
popu lation (|Bergvall fc 6stlidl20ol IZackrisson et al.ll2006l . 
2009). The measured colors could be modeled with a low- 
metallicity (Z = 0.001) stellar population with a very 
bottom-heavy IMF ( dN/dM oc M~ a with a = 4.50) and 
some dust reddening. IZackrisson et all (|2006t ) found that in 
addition to these BCGs, the red col ors of the fai n t out - 
skirts of the stacked edge-on disks of IZibetti et all (|2004T ) 
could also be explained by such an IMF. This soluti on also 
works for the lBergvall. Zackrisson fc Caldwelll (|2009l ) detec- 
tion, sin ce the halo colors of their SDSS stack are similar to 
those of IZibetti et all l|2004l ). If such red halos are found to 
be a common feature of galaxies of very different morpholog- 
ical types, this could substantially contribute to the bary- 
onic dark matter content of the Universe, supporting the 
idea that a non-negligible fraction of the missing baryons is 
locked up in faint low-mass stars. These low-mass stars can- 
not be smoothly distributed over the halo, but would instead 
need to be concentrated in star clusters, since a smooth halo 
component with a bo ttom-heavy IMF is ruled o ut by Milky 
Way star count data (IZackrisson fc Flynnll2008l ). 
After the improved implementation of the th ermally puls- 
ing A GB phase in stellar evolutionary models (|Marigo et al.l 
bOQct ) all BCG galaxies in the sample, except Haro 11, could 
be reconciled with a stellar population with a Salpeter IMF 
of an intermediate metallicity of Z = 0.0 04-0.008, or 20- 
40% of Z Q ~ 0.02 (|Zackrisson et alj|2009h . with the caveat 
that these inferred stellar metallicities of the LSB hosts 
are higher than the observed gas metallicities of the cen- 
tral starburst regions (BO02, Z = 0.002, 10% of Z Q ). Haro 
11 had the reddest optical/near-IR color from the sample, 
V — K = 4.2 ± 0.8, which could not be expla i ned w ith a 
Salpeter IMF. Figure [7] shows the lMarigo et alj |2008h stel- 
lar evolutionary tracks for all galaxies in the sample for a 
Salpeter IMF. Haro 11 is clearly an outlier and cannot be 
reconciled with any reasonable metallicity. The only model 
that provides a reasonable fit for V — K — 4.2 ± 0.8 is the 
botto m-heavy IMF model, al beit with somewhat high metal- 
licity (jZackrisson et alj |2006). 

Analysis of the deep new data presented in this paper sug- 
gests that the color of the host galaxy of Haro 11 is instead 
V-K = 2.3 ± 0.2. This re vised color is shown in Fig. [8] 
with the lMarigo et all l|2008l ) stellar evolutionary tracks for 
a Salpeter IMF for low and intermediate metallicity. We do 
not have new B band data for Haro 11 which would allow 
us to reduce the uncertainty in the B — V color, therefore 
we cannot give a better age estimate for the LSB host. How- 
ever, we have reduced the uncertainty in the V — K color 
by a factor of 4 and from Fig. [H] we can draw the conclusion 
that the host galaxy is consistent with a metal-poor, very 
old, but normal IMF stellar population. 



4.1 Optical/near-IR color interpretation 



4.2 Comparison with previous observations 



The previous color measurement of the LSB host of Haro 
11, V — K = 4.2 ± 0.8, was part of a larger deep 
optical/near-I R photometric study of four luminous metal- 
poor BCGs jBergvall fc Ostlinl 120021) and an additional 
six metal-poor BCGs (|Bergvall et al l l2005h . The colors 
for all galaxies in the sample were found to be consis- 
tently too red to reconcile with a normal stellar popula- 



The discrepancy between our and the BO02 color mea- 
surement is ~ 2cr. The source of this discrepancy is most 
likely the improved quality of CCDs and NIR arrays, as 
well as the larger detector sizes of the new observations. 
The most dramatic change in data quality is seen in the K 
band (Figs. [TJ [2]). Because of the higher sky brightness and 
the lower atmospheric transmission in the K band, ground- 
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based V band observations are usually ~ 5 (Vega) magni- 
tudes deeper than ground-based K band data. The behavior 
of a V — K color profile at very faint isophotes is thus al- 
ways subject to limitations imposed by the depth of the K 
band data. Given favorable observing conditions and assum- 
ing the reduction is carried out correctly, this in turn would 
depend on the detector quality, the total exposure time, and 
the specific morphology of the target galaxy. Since the mor- 
phology of Haro If has not changed in the last decade, we 
examine the effects of the increase in exposure time. For the 
K band data the effective exposure time (taking telescope 
aperture into account) is now 2.4 times longer than the pre- 
viously published data. Comparison of the two K surface 
brightness profiles in Fig. [5] shows that the added detector 
integration has increased the radial range within which the 
K profile is accurate to more than ±0.5 mag by ~ 2 kpc. 
To test whether the longer exposure time also accounts for 
the different slope of the profiles at large radii, we compared 
a K band surface brightness profile obtained with f 00% of 
all new K band frames to a profile obtained with only 40% 
of the frames, equivalent to an effective exposure time of 
25 minutes (25 minutes on the ESO NTT correspond to 
66 minutes on the ESO 2.2m). We conclude that there is 
no indication that a shallower profile would be systemat- 
ically brighter at smaller radii than a deeper profile. The 
effect is simply an increase in the random noise at smaller 
radii but the overall shape and slope of the profile at faint 
isophotes show no significant change. On the other hand, 
an insufficient sky subtraction would cause a brightening in 
the surface brightness profile at large radii from the center 
of the galaxy. To test this, we have simulated a sky sub- 
traction problem by adding a constant sky residual to each 
data point of the ESO NTT SOF1 profile seen in the right 
panel of Fig. [5] We find that in the framework of the ESO 
2.2m IRAC2 data, a sky residual on the order of ~ lAo is 
needed to account for the observed brightening of the ESO 
2.2m IRAC2 surface brightness profile seen in the left panel 
of Fig. [5] We therefore conclude that the shape of the ESO 
2.2m fRAC2 K band profile is most likely due to the lin- 
gering presence of a sky residual of ~ f Act. Over the same 
radial range, this sky subtraction problem has been avoided 
in the ESO NTT SOFI profile mainly due to the improved 
observing techniques, the better quality of the NIR array, 
improved sky region estimation techniques and the use of 
a larger light-gathering area. Specifically, the effective field 
of view of the final stacked image has increased by a factor 
of 2.3, thus providing a comparatively larger area for sky 
estimation, well away from the faintest outskirts of Haro f 1 
visible in the K band. 

The work in this paper indicates the need for a possible re- 
vision of t he color measu r emen ts for the fO BCGs in the 
BO02 and iBergvall et ail (|2005l ) sample. Currently under- 
way is the reduction of a sample of nearly 40 BCGs in the 
optical and NIR, which also includes new and deeper data 
for 5 of the 6 BCGs in the BO02 sample. 



5 SUMMARY 

We have presented deep new optical/near-IR data of the 
BCG Haro ff and performed surface photometry on the 
faint outskirts of the underlying LSB host galaxy. We have 



obtained a total color of V — K — 2.3 ± 0.2 for the host 
over the range r = 5-fO kpc. This is 2.3a off from the pre- 
viously published measurement of V — K — 4.2 ±0.8, which 
was measured out to r = f2.3 kpc. Increasing the radial 
range of the total color measurement out to r = 12.3 kpc 
gives V — K — 2.4 ± 0.5 which has a ~ 2a discrepancy to 
previous results. The new color places Haro II on the blue 
end of the V — K range occupied by the other BCGs in 
the sample and reconciles the LSB host of Haro II with a 
metal-poor (Z = O.OOf) stellar population with a standard 
Salpeter IMF. Hence, the metallicity indicated by the color 
is consistent with that measured from emission-line ratios, 
and an anomalous stellar population is no longer required 
to explain the properties of the Haro II host galaxy. 
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